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' We study analytically the asymptotic behaviour of the average probability V(n,t) for the tra- 

i jectory of a 2D Brownian particle wandering in the presence of randomly distributed traps to wind 

■ n times around a given point after a time t. It is shown that V(n, t) ~ exp(— c-\/t)(l + x 2 ) -1 with 

— . ' x ~ njyi, where the first exponent represents a well-known long-time tail of the probability that a 

particle will not be trapped. 

m ; 

The properties of random walks with various topological constraints has attracted a great deal of theoretical 
interest for many years. Apart from apparent practical relevance to the physics of polymers or Abrikosov vortex 
lines in superconductors (see, e.g., Ref. M), studies of such random walks make profound and intimate connections 
• to many beautiful mathematical results. Perhaps, the most prominent example is the problem of the winding angle 
| distribution in two dimensions. The winding angle of a planar random walk is, by definition, the total continuous 

S angle 8(t) — 2nn(t) swept by a Brownian particle around a prescribed point after a time t. It was found by Spitzer 
H that the asymptotic probability to wind n times is given by a Cauchy law: 
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S: p (njt) ^ 1 s „Jt, a t^oo. (1) 
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d . This result was later confirmed by many authors by employing different techniques (see, e.g., Refs. |3|-|5[]). In general, 
one could also ask how many times a particle has wound around a set of N prescribed points || or an excluded disk 
in a 2D plane §• 



Let us now suppose that our Brownian particle can not wander freely, but instead can be irreversibly trapped 
by the impurities located at some randomly distributed points in a plane. It is known that the properties of such 
a system differ drastically from those of an ideal random walk. For instance, the survival probability is given by 
P{t) ~ exp(-c v / t) §, while the mean square displacement is sub-diffusional: (r 2 ) ~ y/t 0. The purpose of the 
present article is to discover how the presence of traps affects the winding angle distribution function. 

The probability distribution for a random walk starting at a point r' to end at a point r after a time t satisfies the 
diffusion equation 



BP 

^^DV 2 P~U{y)P. (2) 

Here D is the diffusion coefficient and U(r) = Uq J2i <K r — i s the random "potential", which is the probability per 
unit time for a particle to be trapped (Uq > 0). The positions of point-like traps are distributed uniformly in a 
plane according to the Poisson law with mean density p. The solution of Eq. (g) is given by the Wiener path integral 
' formula: 



P(r, t; r', 0\U) = J* _J Vv(t) exp j- J dr (^r 2 (r) + U(v(t))\ j . (3) 

Ch ! r(0) ~ r 

The probability for a closed random walk with r = r' of "length" t to wind n times around the origin in a given 
distribution of traps can be calculated by inserting a (5-function constraint || in (^|), so that 



X • V(n,t\U) = (5(n-±- / dr 9{r)\) , (4) 

f-H \ \ Z7T JO / / P(r,t;r,0\U) 

' 

where 8(t) is the angle between the radius-vector r(t) and some fixed direction in the plane (note that n can be 
non-integer). Writing the (^-function as an integral over an auxiliary variable p, we arrive at 

V(n,t\U)=J dpe 2 * ipn £ {t Pr(r) expj-^ dr ^r 2 (T) + C/(r(r))+i^(r)j|. (5) 

If one assumes that the starting point r is not fixed, then the path integral on the right-hand side is nothing but 
the partition function at inverse "temperature" 1/T = t of a particle of unit charge and mass m = (2D)^ 1 moving 
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in a random potential and in a solenoid field localised at the origin, the solenoid carrying a flux (f> = —2np. For the 
average probability we then have 

P(n, t) ee (V(n, t\U))u = j°° g J™ dE e~^ n e- Et N{E, 0), (6) 

where N(E, </>) is the average density of states. Since we are interested in the asymptotic behaviour of V{n 1 1) at large 
t, all we have to do is to calculate the asymptotics of N(E, 4>) at small E, which is called "Lifshitz tail" flCf |. 

From the analysis of a random walk with traps but without solenoid it is known that at E — ► the main contribution 
to the density of states comes from the large regions in real space which are free of traps. The probability to find such 
a region of area S is exponentially small: p(S) ~ e~ pS . From the elementary quantum mechanics we know that the 
ground state energy of a particle in a 2D potential well with radius R is given by E(R) ~ DR~ 2 ~ DS^ 1 . Therefore, 
S(E) ~ D/E, and the density of states is N(E) ~ p(S(E)) ~ exp(— const pD/E). More formally, such exponentially 
small tails of the density of states correspond to the contribution of instantons jy], which are spatially localised 
solutions of the saddle-point equations in the functional-integral representation of the problem JlJ,[l3| . Our strategy 
is as follows. First, we formulate the problem in the language of quantum field theory with some effective action S. 
Then, we find an explicit form of the instanton solution. The last step is to calculate the asymptotic behaviour of 
V(n,t) at large t due to the instanton contributions. 

The Schrodinger equation for a quantum particle moving in the field of a solenoid and in the potential U(r) is as 
follows: 

Hip ee Z)(— iV — A(r)) 2 '0 + U(r)ip — Eip, (7) 

where Ag = cj)/2irr is the vector potential created by the solenoid. The density of states is proportional to the 
imaginary part of the Green function Ge(j, r) = (r\(E — H + i0) _1 |r), which can be calculated by standard means of 
the quantum field theory. Using the replica trick, the disorder average of the (Euclidean) generating functional can 
be performed, and we have in the limit n — > 

G B (r,r) = J PVW e- 5[V(r V(r)^i(r), (8) 

where <p is an n-component Bose field and V 2 (p = n"=i D(p a 'D<p a . The action is 
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J d 2 r {<p(-E + D(-iV - A) 2 )^ + p(l -e~ Uo ^y (9) 



At E < the action is always positive, so that the field theory is stable and the imaginary part is zero. At 
< E <C E c = pUo (E c being the mean value of the random potential) there is a metastable vacuum state ip = 0. In 
this case, the small- E asymptotics of the imaginary part of the Green function is determined by a non-trivial saddle 
point of the action fPH] and, with exponential accuracy, 

Due to the rotational symmetry of Eq. (|^) in the rt-dimensional replica space, the saddle-point solution (instanton) 
has the form 

<Pa(r) = <^(r)e a , (11) 



where e a is the a-th component of an arbitrary n-component unit vector. From (|9|) and (|l_y) we obtain the following 
equation for the function (p(r) which is assumed to be rotationally invariant in real space (i.e. y(r) = (p(r)e 1 with 
to = 0): 

-D~(r4-)'P + D^ + E c e- u ^' t p = Eip, (12) 
r dr \ ar J r z 

where v = \<p\/2iv. Let us introduce the dimensionless variables: 

r = &, tp(x) = U ~ 1/2 f{x). 

Here £ 2 = D/E is a characteristic scale of the problem, which is nothing but the typical length of diffusion in time 
t = E -1 . Eq. (|l2|) can then be written as 
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--T \ x l~ / + + - / / = /, 13 

x dx \ dx J x^ 

where a 2 = E c /E. 

Since Eq. ( |l3| ) is a non-linear differential equation, we are able to find only an approximate solution. As shown in 
Appendix, at E <C E c (i.e. a ^S> 1) one can replace the "potential" V(f) = a 2 e~f in Eq. (13) by the potential well 
with infinitely high walls, having the shape of a coaxial ring with the inner and outer radii x\ and X2 respectively. 
Then the instanton solution inside the ring satisfies the Schrodinger equation for a particle in the solenoid field, so 
that 

f(x) =A 1 J u (x)+A 2 Y v (x), &tx!<x<x 2 , (14) 

where J„ (x) and Y v (x) are the Bessel functions of the first and second kind respectively. The positions of the matching 
points Xi(y) and x 2 (v) are to be determined from the following equations: 

XiF u (xi)J v {xi) = -X 2 J v {x 2 ), 

xiF v (xi)Y u (xi) = -x 2 Y u (x 2 ), 
where 

I' v {axi) v I v+ x{axi) 



F v (x x ) 



I v (axi) axi I v (axi) 



(the details of derivation can be found in Appendix). These equations are valid for E <C E c and arbitrary v. Going 
back to the dimensional variables, it is easy to convince oneself that the instanton action coincides with the area of 
the ring: 

S^t = 2tT(0 / dr r = it p£ 2 (x\(v) - x\(v)). (16) 

Jri 

In the absence of solenoid (i.e. at v = 0), the solution of Eqs. (|l5| ) is x\ — 0, x 2 — a, where a w 2.405 is the first zero 
of the function Jo (a;). After substitution in (|l6|), the Lifshitz result [^| is recovered. 

At v 7^ 0, due to fast oscillations of the phase factor er^ n in (||), the main contribution to the integral comes from 
small 4>, which allows one to use a perturbative expansion in powers of v. We seek a solution of (jll^) in the form 

xi=5xi(v), x 2 = a + Sx 2 {v) [ol6x\ — * 0). (17) 

The Bessel functions can be expanded in powers of their index Jlifl : 

rjr TJ 7TZ/ 

J v {x) = J (x) + — Y (x) + 0{v 2 ), Y v (x) = Y (x) - —J {x) + Oiy 2 ). 

However, one should be careful in dealing with such expressions, since the Bessel functions are not analytical at x — 0, 
so that we are able to safely expand only the right-hand sides of Eqs. (|l5|). Using the small- a; expansions of the Bessel 
functions on the left-hand sides, we obtain, in the leading order in 5x\ y2 and v. 

1 1 ( aSxtY ira 

aJi{a)dx 2 — y ( a )^i 



aT(u) V 2 J w 2 

1 T(u+ 1) / a8xi\^ v 



a it \ 2 

where T(x) is the Gamma function. Therefore, the solution of Eqs. (|15[) looks as follows: 



2 {TraY Q (a)a\ 1/u _ 2 



a V r(i/ + l) 7 a ( 18 ) 

X2 a+ \ 2Jria) + 7Ta 2 J 1 (a)Yo{a)a 2 ) 1/ ~ a + 2J x {a) 

At ^ — > xi vanishes faster than foa, so that its contribution to the instanton action can be neglected. From (|l6|), 
we then obtain at 6 — * 0: 
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S inst (E, 0) = (1 + b\</>\ + O(0 2 )) , (19) 

where b = 2^!°!) ~ 0-204. A non-analytical dependence on the magnetic flux is related to the fact that one can not 
regard the solenoid field as a small perturbation due to r~ 2 -singularity at small r. 

Finally, we see from ([h]) that the asymtotics of the average density of states in the presence of solenoid is given by 

/ irpDa 2 (l + b\6\)\ , , 

N(E, (/>) ~ exp (—£ K - . (20) 

To calculate the pre-exponential factor, one should make an expansion around the instanton configuration and integrate 
over all non-zero modes. We, however, shall not proceed in this way further and restrict ourselves by the exponential 
accuracy. After substitution of ( p0| ) in (||), we arrive at 

7>(M)~ f# f^e^e--exp^^i±M 



E 

At large t the integral over E can be calculated by the method of steepest descent, resulting in 



P(n,t)~ f d(f> e- i4 > n e^(-2a^TT P Dty/l + 

J — OO 

v -1 



~ exp(-2 av ^Dl) \l + —) , (21) 

where c = ira 2 b 2 f=a 0.756. In calculating the last integral we used the fact that the main contribution comes from 
small values of the flux, which justifies an expansion in powers of <p. The exponential factor on the right-hand side of 
( pl| ) represents the asymptotic probability for a particle without solenoid to survive after a time t and coincides with 
the result of Balagurov and Vaks ||] . The second factor can thus be interpreted as the conditional probability for a 
particle which has survived to wind n times around the origin. 

It is expedient to compare our results with what is known for other similar systems. For an ideal random walk 
without traps, the scaling variable is x = n/lnt, whose asymptotic distribution is given by Spitzer's law (if). For 
a self-avoiding random walk without traps, the scaling variable x = n/ \/lnt has a Gaussian distribution |15,|l6|]. In 
that case, due to the hard-core repulsion, the trajectory wanders farther away from the origin than does an ideal 
random walk, which reduces the winding number. In our case, we see that the scaling variable is x = n/yt, and the 
asymptotic distribution obeys a Cauchy law. The increase of the winding number can be qualitatively understood as 
follows. We are considering the conditional probability, which implies that the particle has survived until the time 
t. This, in turn, means that it spent much of its life in a finite region of the plane almost free of traps and thus has 
never wandered too far away from the starting point. Such a restriction obviously results in increasing entanglement. 

The author would like to thank D. E. Khmel'nitskii for numerous stimulating discussions. The work was financially 
supported by the EPSRC grant no RG 22473. 



APPENDIX: 

We replace the "potential" V(f) = a 2 e~f 2 in the nonlinear equation ( p^|) by a piecewise constant effective potential: 

v eS (f) = {f ;^<J; (ad 



Then the solution is a piecewise continuous function: 



fi(x) , at < x < x\ 1 
f( x ) = { hi x ) > at xi < x < x 2 , (A2) 
f 3 (x) , at x 2 < x, 



where the functions fi(x) obey the following linear equations: 
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I /l,3 + — /l,3 + a /l,3 = /l,3, 



X" 
,2 



x da; 

-IjL (x—\ / 2 + — / 

x dx \ dx 1 x 2 



(A3) 



/a 



The solution /(x) and its derivatives must be continuous functions of x at x = £1.2 ■ The positions of the matching 
points X\ and xi are de term ined from the conditions fi(xi) = f 2 (xi) — 1 and f 2 (x 2 ) — fs(x2) — 1 (see Fig. |l|). 
The solution of Eqs. (A3) looks as follows: 



= C\I v (ax), 

--Ci 1} J„{x) + ci 2) Y„(x), 
h = C 3 K u (ax), 

where I v (x) and K v {x) are the Bessel functions of imaginary argument. The boundary conditions read 

h(xi) = h(x\) = h h{x 2 ) = h{x2) = 1, 



(A4) 



(A5) 



I 1 \ K (1 2) 

Substituting (A4) in (A5), we obtain a system of six transcendent equations to determine C\, , C3, X\ and x 2 - 



fi 2) 

After changing notations C\ — aAi t 2, the equations for Ai^ and X1.2 take the form 



Ai J v {xx) + A 2 Y v {xi) = a- 1 , 
A!J u {x2) + A 2 Y v (x 2 ) = a' 1 , 
A 1 J'M+A 2 Yl(x 1 ) = ij^, 
A 1 J' v {x 2 )+A 2 Yl{x 2 ) = §^. 



(A6) 



In the limit a ^> 1 the right-hand sides of the first two equations vanish. If one assumes that x 2 ~ 1, then in the 
same limit the right-hand side of the last equation tends to —1. Excluding A1.2 from (|A6|), we arrive at Eqs. @. 

It is also worth explaining why we choose to split the plane into three different regions. If there were no solenoid 
(y = 0), then x\ — and the instanton solution would be given by the Bessel function Jq(x) which tends to a constant 
at x — * 0. However, if v 7^ 0, then the naive assumption that x\ = and f(x) ~ J v {x) is not consistent with the 
condition that f(x) > 1 everywhere inside the potential well, since J v {x) ~ x" at x 0. For this reason one has to 
introduce the inner matching point x\ 7^ 0. 
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